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Abstract 

Let X be a linear diffusion and / a non-negative, Borel measurable 
function. We are interested in finding conditions on X and / which 
imply that the perpetual integral functional 



is identical in law with the first hitting time of a point for some other 
diffusion. This phenomenon may often be explained using random 
time change. Because of some potential applications in mathematical 
finance, we are considering mainly the case when X is a Brownian 
motion with drift > 0, denoted {B'f^ : t > 0}, but it is obvious 
that the method presented is more general. We also review the known 
examples and give new ones. In particular, results concerning one- 
sided functionals 



are presented. 

This approach generalizes the proof, based on the random time 
change techniques, of the fact that the Dufresne functional (this cor- 
responds to f{x) = exp(— 2x)), playing quite an important role in the 





and 




study of geometric Brownian motion, is identical in law with the first 
hitting time for a Bessel process. Another functional arising naturally 
in this context is 

roo 

/ (a + exp(4^)))-'dt, 
Jo 

which is seen, in the case fi = 1/2, to be identical in law with the first 
hitting time for a Brownian motion with drift = a/2. 

The paper is concluded by discussing how the Feynman-Kac for- 
mula can be used to find the distribution of a perpetual integral func- 
tional. 

Keywords: Time change, Lamperti transformation, Bessel processes, 
Ray-Knight theorems, Feynman-Kac formula. 

AMS Classification: 60J65, 60J60, 60J70. 



1 Introduction and summary of the results 

Let B^f^^ = {B^^^ := Bt + fit : t > 0} be a Brownian motion with drift 
/i > and / a non-negative measurable function. Encouraged by a number 
of examples listed below, we wish to gain better understanding when an 
integral functional of the type 

POO 

Uf):= / fiBi^^)ds 
Jo 

is identical in law with the first hitting time of a point for some other diffusion. 
Clearly, we can pose an analogous question for an arbitrary diffusion instead 
of B^'^\ In fact, the results in Section [2l are fairly easily extended for arbitrary 
transient diffusions determined by a stochastic differential equation. Our 
interest in the particular case with B^^^ is motivated by the numerous studies 
and results associated to the functional 

POO 

/ exp {-2Bi>'^) ds. (1) 
Jo 

This functional was first considered by Dufresne in [12] where it is seen, 
among other things, how the functional arises as a perpetuity after a 
limiting procedure in a discrete model. 
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We now review some cases of perpetual integral functionals which are 
identical in law with the first hitting time. Let 

Ha{Z):=ini{t: = a] 

denote the first hitting time of the point a for a diffusion Z. 

1) In Yor (see IHj for an English translation) it is shown that for the 
Dufresne functional ^ we have 

r exp{-2Bi>^^) ds Ho{R^'^), (2) 

Jo 

where i?*^^^ is a Bessel process of dimension 5 = 2(1 — /i) started at 1, and — 
reads "is identical in law with". Recall also that 

reM-^Bi^^)ds (3) 
Jo 

where 7^ is a gamma-distributed random variable with parameter /i. We 
refer to Szabados and Szekely for a discussion of Dufresne's functional 
for random walks. 

2) The Ciesielski-Taylor identity: 

j\^^,s^2,^^^ds H,iR^'^), (4) 

where R^^'^ and R^^^"^^ are Bessel processes of dimension 5 > and 5 + 2, 
respectively, started at 0. For a proof, see Williams |J2j p. 159 and 211, and 
Yor (in], im p. 50; in the case 5 = 1 there is a pathwise explanation due to 
D. Williams. We refer also to Getoor and Sharpe p. 98, and to Biane 
for a generalization to a vast class of pairs of diffusions. 

3) The identity due to Biane [2j and Imhof |T9] : 

r^Bi^Uoyd^ - MB^"^), (5) 



where A is a random variable independent of B^^^ and exponentially dis- 
tributed with parameter 2/i and B^ = 0. 

4) The identity due to Donati-Martin and Yor ^ p. 1044: 

Jo exp[2Rs ) — 1 
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where R^^^ is a three-dimensional Bessel process started from 0. See [E] also 
for a probabilistic explanation of 

In Section [21 of this paper we present a general method based on Ito's and 
Tanaka's formulae and random time change techniques which connects the 
distribution of a perpetual integral functional to the first hitting time. This 
method gives us the identity ^ and also the following (which could be called 
the refiecting counterpart of Q): 

^ exp(-2a5(^))l^^(,)^Qjrfs H,/a{R^'^), 6 = 2fi/a. (7) 

where the Bessel process i?^^^ is started at and, in the case < 5 < 2, 
refiected at 0. However, the simplest case emerging from our approach leads 
us to the identity 

r{a + exp{Bi'/'y)r'ds = Hr{B^^'^^), (8) 
Jo 

where a > 0, B^^"^^ = and r = -log(l + a). The reflecting counterpart of 
© is 

Jo [a + exp(iJs )) 

where a > 0, B^"-/'^'> is reflecting Brownian motion with drift a/2 started from 
and r is as above. In Section 3.3, when analyzing the functional on the left 
hand side of (jHl), we also flnd the diffusion with the flrst hitting time identical 
in law with the functional 



oo 

(a + exp(5(^)))-2rfs. 



where /i > is arbitrary. The Laplace transform of this functional can be 
expressed in terms of Gauss' hypergeometric functions (see |3]). 

We have not investigated or constructed a discrete model (as is done in |T2] 
for the functional in (jlj) which would lead to the functional in (jH]). Notice, 
however, that 

/"OO />oo 

/ (l + exp(5(^)))-2rfs = / exp(-2 5(^))(l + exp(-5('^)))-2cis 
Jo Jo 
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and, hence, this functional can be considered as a modification of Dufresne's 
functional such that the discounting is bounded (we call this modified func- 
tional a translated Dufresne's functional). In fact, using the results in Salmi- 
nen and Yor ^Hl, where the integrability properties of perpetual functionals 
are discussed, it is seen that while Dufresne's functional does not have mo- 
ments of order m > fi (cf. Q), which is perhaps unrealistic from an eco- 
nomical point of view, the functional in Q has some exponential moments - 
being in this respect more appropriate. 

For functionals restricted to the negative half line we cannot in general 
have similar descriptions in terms of first hitting times. A typical example is 
the identity ^ above. However, the Lamperti transformation allows us to 
connect exponential functionals to the occupation times for Bessel processes. 
In Section im we show the identity 

/•oo poo 

/ exp(-2aB('^)) = / l^^,s,,^^^^^ds, (10) 

«/ •/ 

where R^^^^ is a Bessel process of dimension 62 = 2(1 — ^) started at 1/a. 
Notice that i?*^^^) j^j^g g in finite time, and, in case < /x < a, we take to 
be a killing boundary point. Further, also by the Lamperti time change, we 
have 

/•oo poo 

/ exp(2aS('^)) = / l^^^s.^^^^^^ds, (11) 

•-'0 

where R^^^^ is a Bessel process of dimension ^3 = 2(1 -|--) started at 1/a. This 
identity was first observed and proved in Yor jig] (^3 P- ^33) by different 
methods. 

Finally, we recall the recent works by the second author, jointly with H. 
Matsumoto, see [21], IHOI, IHII, in which the variable 

[ exp(-25(^)) ds, fi>0, (12) 
Jo 

plays an essential role in obtaining a variant of Pitman's theorem. Indeed, it 
is proved that the process 

Zi'''> := exp(5;^^) [ exp(-25(^)) ds, t > 0, 
Jo 

is a time homogeneous diffusion. Other properties of the functional in (fT2|l 
are studied in the papers [2H1, (E], [Z|, |H], [H]- See also Dufresne [iH] and 
Matsumoto and Yor [321 • 
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The paper is organized as follows: In the next section a general method 
connecting the distribution of a perpetual integral functional to a hitting 
time is presented. Examples of the method are presented in Section |H1 In 
particular, we discuss the functionals (and their reflected counterparts) ap- 
pearing m Q and ©. Further, using the Lamperti transformation we derive 
the identity (jH) from ^ and give a new derivation (for another one, see (SHI) 
for the joint Laplace transform of the functionals 



where L^(i?(^/^)) is the ultimate value of local time at oi B^^^'^\ Especially, 
this latter computation has interesting connections to some earlier works. We 
also prove in Section 01 the identities (|i6|) and ()T7|l (see the table below). 
In Section |3l we modify the Feynman-Kac formula to be directly applicable 
for computing the Laplace transform of a perpetual functional and discuss a 
characterization due to Biane [5] for one-sided functionals. We flnish with a 
short Appendix containing the Ray-Knight theorems used in the paper. 

To summarize the discussion made throughout this work, we use system- 
atically random time changes in a set-up which a priori encompasses the 
scale and speed Feller type representations of one-dimensional diffusions (for 
which, see, e.g. the recent paper by McKean ISBl) as well as Lamperti's trans- 
formation and Ray-Knight theorems. Of course, this stochastic approach 
and the results it allows to derive agree with the more analytic Feynman- 
Kac approach of solving ODE's, and performing for them the corresponding 
(deterministic) changes of variables. 

We conclude this introduction with a table containing (most of) the func- 
tionals and the associated hitting times discussed in this paper. We use the 
notation B^^^^ for a Brownian motion with drift /i > 0, R^^^ for a Bessel 
process of dimension S, and B^^^ for a reflected Brownian motion with drift 
/i > 0. All the processes in the table, if nothing else is said, are started from 0. 
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Ref. 


Functional (a > 0, /i > 0) 


Hitting/occupation time 


m 






m 


exp(-2a fii^^) c/s 




m 


jg exp(-2a5r ) l^^M^Q^ ds 


Jo l{Rf-2'^/'"'>l/a} 


(ED 


fOO ( I ^ 

Jq exp(2a5rO 1|5(m)<o| ds 


/o" l{4'+"'/'"'<l/a} 

^(2+2^/a) _ , 


(EI) 


r°°(a + exp(Ei^/^^))-2(is 


J. J. J. yj^ J , 

r = ^ log(l + a) 






Hx{B^^^), X ~ Exp(2/i) 




/o°°l{i?(*+2)<i} 




® 


/J"(exp(2i?i'))-l)-icis 




(BED 


/o°°exp(-2i?fV5 




(BID 


Ji^il + R^J'Y'^ds, 7>1 


i7l/(,„l)(i?«27-l)/(7-l))) 
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2 Perpetual integral functionals as first hitting 
times 

Let {B^^^ : t > 0} be a Brownian motion with drift /i > and / a non- 
negative locally integrable function. Then (see Engelbert and Senf |TB| and 
Salminen and Yor jUl] ) 



/^(/) := / f{B[^^)du<oo a.s. ^ / f{x)dx <oo. (13) 



In our first proposition it is stated, under some additional assumptions on 
/, that there exists a diffusion constructed in the same probability space as 
B^^^^ such that the perpetual integral functional /oo(/) is a.s. equal to the 
first hitting time of a point for this diffusion. 

Proposition 2.1 Let / : R i— R be a monotone, twice continuously difFer- 
entiable function such that r := lim^_,oo /(a;) exists. Introduce the additive 
functional 






Assume further that f'{x) ^ and 




Let Z be a diffusion given by 



for t > such that 



at := inf{s : Ig > t} < oo 



Then a.s. 




(14) 



Moreover, Z is a solution of the SDE 

dZt = d(3t + G o f-\Zt) dt, Zo 



/(O), 



where (3 is a Brownian motion and 




(15) 
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Proof To fix ideas, take / to be increasing. By Ito's formula 

= r f{B^s'^)dBs+ l\f\B'f^)fG{Bi^^)ds. 
Jo Jo 

Replacing u by at we obtain 

pat pat 

Zt-Zo^ f{B(/^)dB,+ {f{Bi^^)fG{B(/^)ds. 
Jo Jo 

K={f{Bi^y)f and a[ = ^^{f{Bi^J))-' 

at 

it follows from Levy's theorem that 

pat 

Pf.^ / f{Bi^^)dB,, t>0, 
Jo 

is a Brownian motion, and we obtain the claimed SDE: 

Zt-Zo = (3t+ [\f{Bi^])fG{B^^J)da, 
Jo 

^ Pt+ [ Gof-\Zs)ds. 
Jo 



Because 



Because B^ ' — > oo as i — > oo we have 



lim f{Bl''>) ^r, a.s. 

t—*oo 

It follows now from 

Zit = fiBl'^) < r 

by letting i — > oo that 

-^oo — Hj.[Z^ a.s., 
completing the proof. □ 
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Remark 2.2 Let / and Z be as in Proposition 12.11 Then, modifying the 
proof above, it is seen for x > that 

{r{B^^^)fds = -mi{t: Zt = fix)} a.s. (16) 

We consider next perpetual integral functional restricted on R+. For 
simplicity we take / to be decreasing, and leave the case "/ increasing" to 
the reader. Further, we use the notation 

M^'^^ := sup{5(~^) : s <t} := snp{Bs - fis: s <t}, fi>0. 

Recall that 

{pj^) := Mf"^ - bI'''^ : t>0} (17) 

is identical in law to a reflecting Brownian motion with drift /x, i.e., the 
solution of Skorokhod's reflection equation driven by 

-B^-^^ = {-Bt) + fit, t>0. 

Proposition 2.3 Let f and G be as in Proposition \2.1l Assume moreover 
that f is decreasing and let p^^^ be as given in ifT^) . Let 

At := !\f\pi'^)f ds, and a+ := inf {s : Aj > t}. 
Jo 

Define a process Z via 

Zt:= fiO)-f{p^^}) 
for t > such that < oo. Then 

/+ := rU\Bi^^)f ds inf {t : Z, = r*}, (18) 

^ 

where r* := /(O) — r. The process Z is a solution of the reflected SDE 

dZt = d(3t-Go r\f{0) - Zt) dt + L',{Z), Zo = 0, 

where (3 is a Brownian motion, {L^{Z) : t > 0} is a non-decreasing process 
which increases only on the zero set of Z, and G is as in Proposition \2.1l 
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Remark 2.4 In comparison with the statement in Proposition no- 
tice that the equality in (fH^ holds "only" in distribution. 

Proof The starting point is the observation 

and, hence, we find the law of . As in the proof of Proposition 12. II we use 
Ito's formula 

/(O) - /(p1^^) 

Consequently, for t > such that < oo 



/'(pi'^)) rfi?, - /'(0)m(7) - / (/'(pi^)))^ G(p(^)) 

°« Jo 

/ Go/-i(/(0)-Z,)rfs + L°(Z), 

^0 



where Vl{Z) := -f (0)M^7^ determines a non-decreasing process which 
increases only on the zero set of {Z^ : t > 0}. It holds also that 

0<Zt<r*:= /(O) - r 

for alH > such that < oo. To conclude the proof, observe that 



= lim (/(0)-/(p(^)~ 

°° t— >oo \ 



giving dlHI). □ 



Remark 2.5 The method of random time change is also used by Getoor 
and Sharpe |^ Section 5 as a tool to compute the distributions of stopped 
functionals of Brownian motion. Their idea is similar to the one in Propo- 
sitions 12.11 and 12. 3| that is, to identify the functional with the first hitting 
time of a point for some other diffusion. 
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3 Detailed studies of some perpetual function- 
als 

3.1 Dufresne's functional 

As a very particular case of our method, we study now the identity Q (or 
(fT!H) below) and its various one-sided variants. For more details concerning 
(121), we refer to p. 16. For the joint distribution of the functionals in (0T| 
and (f22|) . see Salminen and Yor [221 • Section 4 Example 4.7 an additional 
characterization is derived for the functionals in (f22|) . 

Proposition 3.1 Let B^^^ with fx > be started from 0. Then for a > 
the following 5 identities hold 



where R^^ ^'"/"^ is a Bessel process of dimension 2 — 2/i/a started at 1/a. 



where jg started at and, in the case < n < a, reflected at 0. 



a: (cf. 0), 




(19) 



b: (cf. 0), 




(20) 



c: (cf. 0) 




(21) 



where R^"^ ^^'/"^ is started at 1/a and killed when it hits 0. 
d: (cf. (fTOl) and ((69|)). 




(22) 



where R^^ 'ii^la) started from 1/a and killed when it hits 0. 
e: (cf. (HID), 




(23) 



where R, 



,(2+2/i/a) 




1/a. 
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Remark 3.2 In fact, as is clear from the proofs below, the identities (fTn|l . 
(|2H), (1221), and hold a.s., i.e., the Bessel processes appearing therein can 
be constructed in the same probability space where B^^^^ is given. 

Proof 

a: We let a = 1; the case with arbitrary a can be treated using the scal- 
ing property (however, see the proof of ((2011 where we work with arbitrary 
a). Apply Proposition 12.11 with f{x) = exp{—x). Now / is decreasing, and 
straightforward computations show that 

G{x) = fi)exp{x). 

The process Z is given by the SDE 

dZt = dPt + ^^^^dt, Zo = /(0) = l. 

Consequently, Z is a Bessel process of dimension 6 = 2 — 2/x, as claimed in 
(CI). 

b: Here we use Proposition 12.31 with /(x) = a^^ exp(— ax) which leads us 
to consider the additive functional 



A+ := [ exp{~2a p^^^) ds 
Jo 



where p^^^ is a reflecting Brownian motion with drift p > (see (jT7|l ). By 
Proposition 12.31 the process Z given by 



Z,:=-(l-exp(-ap|fi)) 



satisfies the SDE 



= A + (- - b f (- - ZsY' ds + L°(Z), Zo = 0. 
a 2 Jq a 

Recall that a"*" is the inverse of A'^, and L^{Z) is a non-decreasing process 
which increases only on the zero set of Z. Proposition 12.31 vields now 



exp{-2aBi^'^)l^^M^^^ds = J exp{-2a pi^^) ds 



13 



It remains to show (cf. (EH)) that Hi/a{Z) is identical in law with 
One way to do this is to compute the Laplace transform of Hi/a{Z) by using 
the standard diffusion theory (see Ito and McKean [20j or 11.10 p. 18). A 
slightly shorter way is to observe that {l/a — Zt : t > 0} is a Bessel process 
of dimension 2 — 2/i/a started at l/a, reflected at l/a and killed when it hits 
0, denoted ^(2-W«). This gives us 

Finally, the claim 

= i7i/,(i?(2M^)), (24) 

where the Bessel process 

j^{2f,/a) -g started at 0, can be verifled by straightfor- 
ward but lenghty computations with Laplace transforms. (As seen in Section 
3.2, the Ciesielski-Taylor identity follows easily from (f2i|) ). 



c & d: To prove (pT|l and (f^ we use the Lamperti transformation (see 
Lamperti jSHI) as in the proof of Dufresne's identity in Yor Notice that 
the Lamperti transformation is used (and proved) implicitly also in the above 
proof of (fT9|) . To recall the transformation, let 



aS^) = f exp{-2a Bi^"^) ds 
Jo 



and 

a[^^ := inf{s : > t}. 
Then the (first) Lamperti transformation says that the process 

|-exp(-a5^';|,) : t>0} (25) 
a "t 

is a Bessel process of dimension 2 — 2/i/a started at l/a and killed when it 
hits 0. Letting /?(2-2/./a) (jenote this Bessel process we obtain 

exp(-2a5^^^) 1|b{m)>o| ds 



oo 



exp(-2a5j^';|)) da'^p 
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oo 







For we have similarly 

exp(-2aS('^)) 1|^(m)<o| ds 

PCX) 

= / exp{-2a B^%) 1 M da'^^'^ 
Jo °= ^ .(1)^"^ 



;.Ho(K(2-2M/a)) 





oo 







e: The identity (f23|) can be proved using the Lamperti transformation ^ 
with negative /i. This can also be formulated for > by defining 



:= /" exp(2a5('^Vs 



and 

a; ' := m 



- ^-f{s : > t}. 



Then the (second) Lamperti transformation states that 

{-exp(a5^'^])) : t > 0} (28) 



is a Bessel process of dimension 2 + 2yu/a, denoted started at 1/a. 

Consequently, as in the proofs of c and d above, 

/•oo POO 

/ exp{2aB^J'^) l^^(^)^^^ds = l^j^{2+2^/a)^^^^^ds. 

«/ •-'0 

□ 
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Remark 3.3 (1) Notice that in case a = 1 and /i = 1/2 the identity 
(see also JSEl)) takes the form 

1" exp(-2 5(1/2)) (A) 

= / l{i?.<i} ds 
Jo 



Ei(exp(^-7^ l{Bs>i}ds^^ = - 



where p is a reflecting Brownian motion started at 0. Similarly, (|77jl can be 
written as 

From PI formula 1.2.4.1 p. 200, for 7 > 

r-Ho(B) 

1;d ^11 dsW = -, 

By the well known formula (cf. [3| formula 1.2.0.1 p. 198) 

E,(exp(-7i/o(5))) = exp(-xv^), x > 0, 7 > 0, 
and, therefore, 

^"exp(-2 5(1/2)) (A) H,iB), 

where B is started at and A is an exponentially (with parameter 1) dis- 
tributed random variable independent of B. 

(2) By the Lamperti time changes (f25|l and (f28|l . a general perpetual inte- 
gral functional of geometric Brownian motion can be expressed in terms of 
Bessel processes as follows: 

/•oo roo 

/ /(exp(aS(^)))cis= / (ai?(2+2A'/'^))-2 /(a^i'''"'/"^) 
Jo Jo 

where i?(2-Wa) jg killed when it hits 0. 
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3.2 Ciesielski-Taylor identity 

The identity (f2n|l (cf. also (f2i|) ) is now applied to deduce the well known 
and puzzling Ciesielski-Taylor identity (jH) (or (f29|) below) , see Ciesielski and 
Taylor 0, and also Yor and gH, Chap. IV. 

Proposition 3.4 The following identity in law holds: 

I l{iji^+^)<i/a}^^ - Hi/a{R^'^), V(5>0 (29) 

where the Bessel processes are started at 0. 
Proof By Proposition 13.11 h and c 

1/a} 

where /^(^"^/i/a) jg g^arted from 1/a. Using the standard time reversal argu- 
ment (see Pitman and Yor ^24j p. 341) we obtain 

/OO POO 
l{iJp-2''/"'<l/a} = l{i?f+2^/'"'<l/a} 

where i?(2+2At/a) jg started at 0. Letting 5 = 2/x/a gives now the claimed 
identity □ 



3.3 Translated Dufresne's functional 

An interesting case emerging from Proposition [O] is when the function G in 
(|T5|l is equal to a constant c, say. Straightforward computations show that 
this yields us the functional 



+ Aexp(2/xS('^)))"^ds, 



where A is a free constant. It is natural to assume that c > and A > 0. By 
the scaling property of Brownian motion, 

Jo Jo fi 
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indicating that we are, in fact, obtaining information only about Brown- 
ian motion with drift n = 1/2. Choosing the constants appropriately, this 
functional, which in [3^1 is called Dufresne's translated perpetuity, can be 
expressed in the form 

POO 

/ {a eMBi^^^^) + l)-^ ds 
Jo 

where a > 0. In [39j we derive the joint Laplace transform for the one sided 
variants of this functional, that is, for the functionals in (f3T| and (f32|l be- 
low. Here we give a new derivation of this Laplace transform based on the 
Ray-Knight theorem for Brownian motion stopped at the first hitting time. 
However, before this, we characterize the functional and its one-sided variants 
via hitting times. 

Proposition 3.5 Suppose that B^^'^^ = 0. Then the following 3 identities 
hold 

a: ^ 

r(aexp(i?(i/2)) + l)-2rf5 (J HripW^)), (30) 
Jo 

where a > 0, r = log((l + a) /a), and P^^^'^^ is a Brownian motion with drift 
1/2 started at 0. 
b: ^ 

{aeMBi'/'^) + iy\^^af^^^,/s H^jS^'l^^ (31) 

where a and r are as above, and is reflecting Brownian motion with 

drift 1/2 started at 0. 

CI 

J^^ (aexp(S(V2)) + iy\^^,,,,,^^^ds H,{P^'/'^), (32) 
where a > —1, A is exponentially distributed with parameter (1 + a), and 

pirn 

is as in 

Proof 

a: Instead of simply refering to Proposition [22 it is perhaps more instruc- 
tive to go through the computation; in this way we also gain better under- 
standing why the case /i = 1/2 is a special one. Define 

f{x) := log (— — — -), 

\a + exp(— x)/ 
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and observe that 

f'{x) = (aexp(a;) + 1)"^ 

By Ito's formula 

log (a + exp(-5j^))) - log(a + 1) = - ^* — 



Introduce 



and 



aexp(5i^^) + l 
/■* ds 

2 Jo aexp(5i^^) + l 2 X (aexp(S^'^^) + 1)^ 

ds 



+ -o / : .„r„.K (33) 



{aexp{Bi^'>) + iy 



Z„ := /(^i'^)) = log ( 



a 



1 

^« = ^0 + + (/i - 2) I Y 



where a is the inverse of A. We have Zq = log(a/(l + a)) < 0, and < for 
all M such that au < 00. Moreover, from (jTH^ it is seen that Z is a solution 
of the SDE 

ds 

- exp(Z^ 

where P^^^'^^ is a Brownian motion with drift 1/2 starting from 0. Conse- 
quently, because ^ as m ^ 00 we obtain, as explained in the proof of 
Proposition 12. H 

dn 

/ = HoiZ), 

Jo (aexp(5i^^) + 1)2 

and, in particular, for = 1/2 we have (f3fl| . 

b: The identity (f3T| can be proved by directly applying Proposition 12 .31 
and we leave the details to the reader. 

c: Consider now the identity (p^ . Our proof of this relies on the Ray- 
Knight Theorem 15.21 given in Appendix at the end of the paper. Firstly, for 
arbitrary /i > we have by the occupation time formula 

(aexp(5(^)) + l)"'l^^(,.)^^^rfs = £ (aexp(y) + i)"'l^(5(^)) dy, 
where L^(i?^^^) is the total local time of B^^^ at level y. By Theorem 15. 2| 



£ {aexp{y) + iy'Ly^{B^^^)dy, (a exp(-|/) + l) Z 



ydy, 
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where Z satisfies the SDE 

dZy = 2^J~Z'ydBy " 2^Zy dl/ . 

The distribution of Zq is the distribution of L'^(B^^^), i.e., the exponential 
distribution with parameter fx. Consider 

(a exp(-y) + Zy-{a + 1)-' Zq 

ry 

= 2 / (a exp(-M) + 1)-^ VZudBu 







f 2^+(2,.-l)aexp(-ii) 



Define 



Cy :-- 



ry 

/ (a exp(— -u) + 1)^^ Zu du, y > 0, 
Jo 

and let c denote the inverse of C. From (|3i|l we obtain for /i = 1/2 

(a exp(-Cj/) + 1)"^ Zcj, - (a + 1)"^ Zq = 2 Py - y, 
where /3 is a Brownian motion. It follows, because — » as t ^ oo, that 

J^^ (aexp(i?(V2)) + = inf{l/ : /^^^/^^ = 

where ,^ = Zo/2{a + 1) is exponentially distributed with parameter (a + 1). 

□ 



Remark 3.6 The formula (fT6|) in Remark 2.2 gives 



/o (aexp(5r/'^) + l)2 ' 

where p > 0, and g = log ^(a + l)/(a + exp(— p)) 



We proceed by computing the joint Laplace transform of the functionals 
appearing in (f3T|l and (|32|l . This Laplace transform is also given in [SHI but 
here we give a new derivation which from our point of view has independent 
interest. Define 

AW:= r ^^"^"J^ ds. 



(aexp(i?, 



(1/2) > 

S I 
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Proposition 3.7 For non-negative k, K and c 

Fik,c,K) := Eo(exp ( - - cLI{B^'^'^) - K Ai"))) (35) 

VSkTl exp(|) 

~ V8k + 1 cosh(| V8k + 1) + (2c(a + 1) + V8K + 1) sinh(| ^/8kTT)' 

where r = log((a + l)/a). In particular, with the notation as in Proposition 

da 

F(A;,0,A;) = Eo(exp(^-A; ^ (aexp(Sf Z^)) + i)-2 ^sj^ 
= Eo(exp(^-fci7,(/5(^/2))^ 

= exp(-^(V8A; + l-l)^ 



F(0,0,ir)=Eo('exp(-K r 

^ ^ Jo (aexp(5f/'^) + l)2 



Eo(exp(-iri7, 
1 + a 



a + I + + i 

and 



F(fc,0,0) = Eo(exp (- A; /" ^ 



(aexp(5F')) + l) 



ds 



\2 



Eo(exp(-fc/7,(/3(^/2)) 



V8A; + 1 exp(^ 



V8A; + lcosh(| + 1) + sinh(| V8A; + 1) 

Moreover, Ai^-* and Ai ■* are conditionally independent given L^{B^^^'^^). 

Proof We begin as in [22] and express the functionals Ai"^-* in terms of a 
Brownian motion without drift. Firstly, notice that 

^(+) (1) /- exp(2i?i-^/^^)l^^(-v.)^^^ 
io (a + exp(5i~'/')))2 
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and 

^(_) (d) /- exp(2i?r^/^^)l^^(-v.)^^^ 

Jo (a + exp(5i"'/')))2 

We apply the Lamperti transformation with /i = —1/2 and a = 1 on the 
left hand sides of these identities. Recalling that R^^'^ is, in fact, a Brownian 
motion killed when it hits we obtain 

where 5' is a Brownian motion started at 1. Instead of working with a 
Brownian motion B' starting at 1 we introduce a "new" Brownian motion 
B = 1 — B' starting at 0. Using spatial homogeneity of Brownian motion we 
obtain 

Jo (a + 5^)^ 

and 

Jo (a + B',y 
By the occupation time formula 

{A\ L^^iB), A-) 

where L^^(S) is the local time of -B at 1 — y up to Hi{B). From the proof 
of (f32|l in Proposition 13.51 we know 

Eo(exp(-irA-) I L%^{B)=u^ 
= Eo(exp(-irA(-)) 

^exp(-^(v/8^-l)). 



(d) 



{Bs>0} 



ds =: A^ 



(a + l-B,Y - • - ' 



(d) /'^^(^) 1 



{B,<0} 

(a + l-B, 



ds =: A" 
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From the Ray-Knight Theorem 15.11 it now follows 
Eo(exp ( -kAM - cL^ (5(^/2)) - K A^^^ 

= Eo(exp( -/c / , \^ ds -cXl^'> ~-^{V8K + 1-1 



where X*^^) denotes the 2-dimensional squared Bessel process started at and 



1 



Next, recall from Pitman and Yor or Revuz and Yor |36| Exercise 1.34 
p. 453 that for a general squared Bessel process X^^\ 5 > 2, and for any 
positive Radon measure m on [0, oo), we have for x > and t > 

E.(exp( - i I^^Xf m(rf.))) =0(tpexpf|0'(O) 



2 



where v ^ (f>{v) is a unique positive, continuous and non-increasing function 
satisfying for < u < v < t 



(f)'{v)-(f)'{u) = (f){s)m{ds), 0(0) = 1, (38) 

J {u,v] 

and which is a constant for v > t. By 0' we mean the right hand side 
derivative of 0. Notice from p8|) that is convex in (0, oo). Choosing 

where A is a Borel set in (0, oo) and £{1} is the Dirac measure at 1, yields 
Eo(exp( - i / Xf)m(ds) 

^ ^ ^ ^(0,1] 

= Eo(exp(-^y^ (a + .)^ 
= 0(1). (39) 



'{0,1] 



ds -7Xf^ 
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To find the function 0, we proceed in a similar manner as in Example 1 
p. 432. It follows from (p?^ that for f < t is the (continuous) solution 
of 

2k 

^"'"' = (^^<^'')- <^°' 

satisfying the conditions (notice that 0'(1) = 0) 

0(0) = 1, 0'(1-) = -270(1). (41) 

It is elementary to check that y{v) = (a + f)" is a solution of pn|l if and only 
if a{a — 1) = 2k, that is 

a = hl± V8k + 1) =: a±. 



Introducing 



y^{v) := {1 + and y_{v) := {1 + -) 
a a 



our task is to find constants A and B such that 

(f){v) := A?/+(i;) + 5y_(t;), 

fullfills (IHJ. We skip the detailed computations and state only the result 
needed in (fHnil : 

0(1) = w fy;(l) + 27y+(l) - ivUl) + 27y-(l))' 



where w = \/8k + l/a is the Wronskian. The desired formula in Propo- 
sition 13.71 results now from and ([HHll by recalling the definition of 7 in 
p7|l and substituting r = log((a + l)/a). □ 



Remark 3.8 It is seen from the proof above (take 7 = 0) that for a > 
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This formula is valid also for < 5 < 2 when the boundary point is taken 
to be reflecting. In particular, X^^^ is a Brownian motion squared. Similarly, 
for a > 1 

Eo ( exp 




-x'4l))-'^\ (42) 
where for v < a 

:= (1 - and := (1 - -)"". 

a a 

The formula (|i2|) is derived in Mansuy |22j for squared Brownian motion 
using different techniques (but it is also indicated therein that the result can 
be obtained in the way presented above). See also Mansuy [211 for closely 
related results. 



3.4 An identity due to Biane and Imhof 

Next we consider the identity ^ (renumbered (|i3|) below) found by Biane 
[2] and Imhof [IH]. This identity is also observed in [SH] in connection with 
a storage process. The distribution of the random variable Hx featuring on 
the right hand side of (HTni is in this context called the RBrownian motion- 
equilibrium-time-to-emptiness distribution, see Abate and Whitt 

We give two proofs of the Biane-Imhof identity (ji3|l . The first one is based 
on the Ray-Knight Theorem 15.21 and the random time change techniques. 
Given these tools the proof itself is very short. In fact, we repeat the idea 
of the proof of (p^ in Proposition 13.51 The second proof is also based on 
random time changes but now we work via Tanaka's formula. This latter 
presentation is close to the one in [TT] (see also Biane |2|), but we wish 
to give it anyway to demonstrate the connections between occupation and 
hitting times. 

Proposition 3.9 For /i > 

\si^)^,yds H,iB^^^), (43) 

where B^^^ is started at and X is an exponentially (with parameter 2^) 
distributed random variable independent of B^^^\ 
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Proof 1 (based on the Ray-Knight Theorem \5.'A ) We use the notation and 
the structure of the proof of (p^ in Proposition 13 .51 Firstly, by the occupa- 
tion time formula and Theorem 15.21 



/*oo /* /* oo 

/ 1 ds= Ly^{B(^^)dy / Zydy, 

Jo J-oo Jo 

where Z satisfies 

dZy = 2 ^/Z^ dBy - 2^iZy dy. (44) 

The distribution of Zq is the distribution of L'^{B^^^), i.e., the exponential 
distribution with parameter /i. For the random time change, introduce 

Cy := I Zudu 
Jo 

and let c denote the inverse of C. It follows from (jiill 

Zcy - Zo = 2l3y-2ij,y, 

where /5 is a Brownian motion. Consequently, because Zy ^ as y ^ oo we 
obtain (|i3|l . □ Proof 2 (based on the Tanaka formula). Consider 

where LP{B^^^) is the local time of B'^^'^ at (with respect to the Lebesgue 
measure). Introduce 

:= j l^^(j,)^^^ds, and a~[ := inf{s : A'^ > t}. 
Then letting i^'"^ = \ Vl{B^^^) 



IjM"' <o} dBs- fx / <o| ds + 



\Bi^Uo}dBs-fxt + &l 
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It is a straightforward application of the result due to Dambis, and Dubins 
and Schwarz (see Revuz and Yor (HS] P- 181) that the process given by 



^' '^ Jo "^{s='''<o}^^' 

is an JF --Brownian motion. Notice that because B^^} < we have for all 
t > 

i^^l>Pt + f^t. (45) 

Clearly, = A^^ for t > Aq := Ao(5('^)) := sup{t : B^^'^ = 0}. Conse- 
quently, {B^^} : t > 0} is defined only for t < A]^^ and 



lim B^''^ = b[^^ = 0. 

^0 



Because a is the inverse of A and i^^^ does not increase after An we obtain 



From (|i5|l it now follows that 

A- = AX„ = inf {t : p, + ^,t = i^^^}. 

Because /? is a Brownian motion and the Po-distribution of is exponential 
with parameter 2fi it remains to prove that f3 and il^ are independent. To 
do this, let 

A^ := J (is, and := inf{s : A^ > t}, 

and proceed as above to obtain 

= ^, + ^t + &}, 
"t 

where {7^ : t > 0} is a Brownian motion. Because 5^+^ is well defined and 
non-negative for all t > we deduce from Skorokhod's reflection equation 
that 

^^''^ - sup {-7s - /is}, 
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and letting t ^ oo gives 

= sup{-7, -/is}. 



s>0 



From the extended form of Knight's theorem (see [SEl P- 183) we know that 
(3 and 7 are independent and, therefore, also and (3 are independent. □ 

3.5 LeGall's identity 

From a representation formula in LeGall [2l| we can deduce the result in 
(ji6|) below, see Donati-Martin and Yor |Tni p. 1044 and 1052. In this section 
we prove (j46|l using time reversal argument and time changes. This nice 
application of the classical time reversal result was pointed out to us by Y. 
Hariya (in the context of (jTTjl ). 

Proposition 3.10 The following formula holds: 

exp ( - 2 i?f ) ) ds = ifi (i? ) , (46) 

Jo 

where the Bessel processes are started from 0. 

Proof For a; > let A^{R^^^) := sup{t : ^ = x}. Then 



/ exp ( - 2 R^) ds = / exp ( - 2 5, 

Jo Jo 



ds. 



where i? is a standard Brownian motion started from x. Let for t > 



At := [ exp {-2Bs) ds 
Jo 



and, as usual, at is its inverse. By the Lamperti transformation ()25)) . the 
process {Zt : t > 0} where Zt := exp(— _Bq,J is a 2-dimensional Bessel 
process. Clearly, 

Zq = exp(— x) and < Zt < 1 W t < Ah^- 

Consequently, 

Ah, = H^{Z), 

and letting a; ^ 00 proves (|i6|) . □ 
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3.6 Hariya's identity 

We learned the identity in the next proposition from Y. Hariya [iHl and offer 
here a proof which differs from Hariya's proof and is in the spirit of the 
present paper. 

Proposition 3.11 The following identity holds: 

ril + R^'Y'-'ds'^ H,/,,_,,iR^% (47) 

where 7 > 1, 5 = (27 — 1)/(7 — 1) and the Bessel processes are started from 0. 

Proof By the scaling property of Bessel processes and the occupation time 
formula the right hand side of (|i7j) can be written as 



(d) 1 



(7-1)^ 

Consequently, by Theorem 15.31 in Appendix (with the notations as therein), 
(|17|l is now equivalent with 



7(2) pi z^'^^ 

^iy'S^S-2)l^iy. (48) 



In particular, when 7 = 2 (and 5 = 3) the identity (ji8|l takes the simple form 
To prove (jlTIl via (j48| we use the well known fact that 



{Zf): s>0}^^ {{l + srz[f^,^^^: s>0}. (50) 

Notice that follows directly from (f5n|) . However, to obtain we have 
to work little more. Substituting x = y^~^ on the right hand side of 
shows that PH| is equivalent with 
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By the representation (fHIHl the left hand side of takes the form 

■7(2) /-oo ^(2) 

^ X (1+1/)^--^ ^' 

and substituting here a; = 1/(1 + ?/) leads to the right hand side of (f5T|l . 
completing the proof. □ 

Remark 3.12 Consider the Ciesielski-Taylor identity (jH): 

/ l{4*+^)<i}^^ - (52) 

where R^^^ is a Bessel process of dimension 5 > started at 0. It is possible 
to express in alternative forms using Hariya's identity and some simple 
transformations. Indeed, by Theorem 15. 21 (b) we rewrite first in the form 
(see Yor [44] p. 52) 

1 7(2) , .1 7(2) 

^„4-2 , (d) 1 r 



6-2 Jo 6 Jo y 

Applying (ji8|l on the left hand side of (f53|l yields 



dy = T / ^^dy. (53) 



(5-2)2 7, iiTw^'y^sl 

Making the change of variables x = y^ on the right hand side of and 
recalling that 7 = (5 — l)/(5 — 2) leads to 



(5 -2)2^0 (l + ^)2(5-l)/(5-2) 52 2.2(5-l)/5 (^5) 

Further, substituting on the right hand side of (f55|l y = 1/x and using the 
time inversion property of Bessel processes we obtain 



1 r°° z 



(2) 

y 



(5-2)2 (1 + ^)2(5-1)7(5-2) 

(d) 1 



dy 



(d) 1 r ' 



^2 ^2(5-l)/5 

) 



(d) 1 r yP 



52 (1 + m)2(5-1)/'5 
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where F^^-' denotes the process Z^"^^ started with an exponential distribution 
with parameter 1/2 (which is the distribution of Zl when started from 0). 

4 Feynman-Kac approach to perpetual integral 
functionals 

In this section we show how the Feynman-Kac formula can be used to find the 
Laplace transform of a perpetual integral functional /oo(/) where / satisfies 
the integrability condition (fT3|l : 

fiy)dy < oo. 

We remark also that in special cases one can find the law of a perpetual func- 
tional by limiting procedures but for a general characterization the problem 
has to be analyzed more carefully. 

For treatments of Feynman-Kac formula, we refer to Durrett and 
Karatzas and Shreve [22]. See also Jeanblanc, Pitman and Yor (SHI for con- 
nections with excursion theory. 

Consider for 7 > (and /i > 0) the function 

X ^ ^,{x) = E,(exp( - /oo(7/))) = E,(exp( - 7^ /(S^'^)) ds^^ 
By the simple Markov property it is seen that the process 

{^,(si^)) exp ( - 7 1 fiBi'^^) dsy.t> 0} 
is a bounded martingale. 

Proposition 4.1 The function is non-decreasing and satisfies for all x 
andt>0 

^,{x) = E.(vI/,(5;^)) exp ( - 7 1 /(5i^)) ds^y (56) 

Moreover, 

lim ^^(x) = 1. (57) 



a;— >oo 
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Proof The formula is immediate from the martingale property. Let 
X < y and apply the optional stopping theorem to obtain 

^,{x) = E.(vI/,(i?Jg)exp(-7^'''/(i?l'^^)rf^)) (58) 

which shows that is non-decreasing. Next notice that 

^.{Bi"^) exp ( - 7 £ rfs) = E(exp ( - 7 /(^l^^) rfs) | -F*) • 

Consequently, by the martingale convergence theorem, 

%{^) exp ( - 7 /(Sl^)) ds) = exp ( - 7 y /(^i^^) ^^s) 
showing ([FTjl . □ 



Remark 4.2 Notice from (j58|] that for x < y 



In many cases the function can be found by solving a second order 
ODE. This is formulated in the following 

Proposition 4.3 Assume that f is piecewise continuous and satisfies the 
integrability condition ifT^) . Then x 1— > "^-yix) is the unique positive, non- 
decreasing and continuously differentiable solution of the problem 

^v"{x) + fiv'{x) --f f{x)v{x) = 0, (59) 
lim v{x) = 1. 

Proof Notice that B^^^ killed according to the additive functional 

/t(7/) = 7 f f{B'f^)ds 
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is a linear diffusion B* (in the sense of Ito and McKean and its basic 

characteristics (speed measure, scale function and killing measure) can be 
determined explicitly (see f3| No. II. 9 p. 17). From (f56|l it follows that 
is an invariant function for B*. It is well known that invariant functions 
of a linear diffusion are continuous (see Dynkin [11] Vol. II p. 7), and 
differentiable when the scale function is differentiable (see Salminen P- 
93). From the representation theory of excessive functions we know that 
for B' there exist two invariant functions, denoted ip^ and ip-y ^ind called 
fundamental solutions or extreme invariant functions, such that if h is an 
arbitrary invariant function then there exists constants Ci > and C2 > 
such that h = ci ip^ + C2 (p-y. Moreover, we have 



P^H^iB') < oo) 



<^^(z)' 



X > z, 
X < z. 



From this representation it follows that ip and ip solve (f59|l on the intervals 
of continuity of /. Next notice that, because /i > 0, 

lim P^{H,{B') < oo) = 

z— >— oo 

giving ^■yiz) +00 as z ^ — oo. Consequently, all invariant non-decreasing 
functions are multiples of ijj^. In particular, is a multiple of and, hence, 
the condition (|H7| determines uniquely. □ 



Remark 4.4 At the points of discontinuity of / the function \Ef^ usually fails 
to be two times differentiable. 

Example 4.5 We compute the Laplace transform of the functional 

/•oo 

/ ia + exp{Bi'/^^))-^ds 
Jo 

appearing in Consider, for a moment, the case with general /i > 0. 
Taking 

f{x) = (a + exp(a;))~^ 
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in (fHn| gives us the equation 

- v"{x) + fiv'{x) — 7 (a + exp{x))^^ v{x) = 0. 
Putting here x = liay and g{y) = v{\ny) yields 

^ y2 g"iy) + (/i + ^) y - 7 (« + Z/)"' giv) = 0, (60) 

which is, of course, the corresponding equation for geometric Brownian mo- 
tion. By Kamke [21] 2.394 p. 497 this equation can be solved for = 1/2 by 
making the substitution 



V{0 = 9iy), e=^ln(^), 
a y + a 

which transforms (f6n|l to the following 

ri" + ar]' = t]. (61) 
Letting P := ajl-J^^ the general solution of (|6T| can be written as 

r/(0 = A exp ( - (VI+/32 + /5)e) + i? exp ((^1 + /^' - /5)e) • 
Consequently, the increasing solution of (|6n|) is 

\ 2p a + y / 

y ^2-1(^/1+^-1) 



^a + y. 

Notice that ^/'(oo) = 1, and it follows 

exp - 7 (a + exp(Sf /2)))-2 

exp(x) \ (2a)-l(i/a2+87-a) 



Ex 



.a + exp(xj 

For a geometric Brownian motion X with Xq = x > defined via 

X. = exp(5(V2)), Bi'/'^=x, 
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(62) 



the formula (fl^ takes the form 

The identity in law in ^ can be deduced from (or Notice also 

that substituting in (fH^ x = and letting a ^ we obtain a special case of 
the identity 



i-e., 

where Z1/2 is a F-distributed r.v. with parameter 1/2. 

Using Propositions 14. Il and l4. 31 we derive an interesting result due to Biane 
[2] which characterizes the law of a perpetual integral functional of B^'^\ jj, > 
0, restricted on R_ in terms of the same but unrestricted functional of another 
diffusion stopped at the first hitting time. We remark that in |^ a more 
general situation (not only B^^^^) is considered. However, the main interest in 
[2] is focused on occupation times, the aim being to generalize the Ciesielski- 
Taylor identity (jij. The result in our Proposition 14.61 is extracted from 
Remarque p. 295 in and formulated for B^f^\ 

Proposition 4.6 Let f be a positive C^-function such that 

f{y)e'^'^y dy = ^, (64) 

and X a diffusion with the generator 

1 f 1 f'{x)\ d 



2dx^ ^ V 2 f{x)) dx 
Then Hq{X) < 00 a.s. if Xq < 0, and, moreover, 

/(i?('^))l^^(.,^„^rf. fiX,)ds 

where Xq is taken to be exponentially distributed on (—00, 0) with parameter 
2fx. 
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Proof Notice that the condition means that the scale function of 
X given by (cf. ^ II.9 p. 17) 

satisfies 

lim S'^{x) = — oo. 

x^—oo 

This implies Hq{X) < +00 a.s. when Xq < 0. From Proposition 14.31 we know 

^,{x) :=E.(exp(-7/oo(/))) 
is the unique, non-decreasing function such that 

=7l(-oo,o)(a;)/(x)^(x) (65) 

and lima;.^+oo ^{x) = 1. For x > we clearly have 

^^(x) = P,(ifo(5('')) = +00) + P,(ifo(5('')) < +00) ^^(0), 

and, hence, it is enough to compute "^-yix) for x < 0. For this, consider the 
equation 

-u"{x) + fiu{x) = ■y f{x)u{x). (66) 

Let V'7 and (p^ denote the fundamental non-decreasing and non-increasing, 
respectively, solutions of and, similarly, 'tp'y and (f^ are the fundamen- 
tal non-decreasing and non-increasing, respectively, solutions of . Notice 
that / does not have to satisfy the integrability condition (fT3|l . However, 
for (j66|) . we can still argue as in the proof of Proposition 14.31 that all in- 
variant non-decreasing functions are multiples of 'tp^. Using continuity and 
differentiability requirements, ip'y can be expressed in terms of 'ip-y as follows 



^Jx) 



ip^{x), X < 0, 
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where S{x) = — exp(— is the scale function of B^^\ Consequently, for 
X < 

and, in particular, for x = 
To proceed, define for y < 

and notice that V^" exists and is continuous because / G Using the fact 
that solves it is straightforward to verify that ip!^ > 0, i.e., is 
increasing, and that ip-y is a solution of the ODE 

^u"{x) + (/i - ^ j^) u'{x) = 7 fix) u{x), x<0. 
By Ito's formula, the process 

rtAHoiX) 

ip^iXtAHoiX)) exp(-7 / f{Xs)ds) : t > 



is a martingale and, further, because it is bounded, we obtain for Xq = x < 
by the dominated convergence theorem 

E.(exp(-7 r^""^ f{X,)ds)) = tM (68) 
^ ^ Jo ^7(0) 

Observe that 

i)^ [x] 2/i e^'^^ dx = 2y. \ (2/iV^^ [x)^^' {x)) e^^"" dx = 2^1^^ (0) . 



Consequently, if Xq is exponentially distributed on (—00, 0) with parameter 
2/i, (inH) and ^ lead to 

/ 2/ie2^^E,(exp(-7 / /(X,) t/sj j t/x 

t/ — 00 J 



Eo(exp(-7^ 



as claimed. □ 
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Example 4.7 Consider the functional 



Recall from (f22|l . that there exists a Bessel process R^^ "^^^ started from 1 
such that 

poo 

1=1 lr„(2-2M)^ii ds a.s. 

As an application of Proposition 14.61 we derive a new charaterization of the 
distribution of /. Taking f{x) = e~'^^ it is seen that the diffusion X in 
Proposition 14.61 is a Brownian motion with drift /i + 1. Consequently, 

I = / exp(-2S(^+^)) rfs, 

Jo 

where B^l'^^^ is exponentially distributed on (— oo, 0) with parameter 2/i. To 
develop this further, let a; > and assume that = -X < 0. We have 



/ exp (-2 ds 

Jo 

/ exp (-2 (5(^+1) - x)) ds 

exp f-2E(^+^M 



e^^ inf{t : i^J-'^^ = e"^}, 

where B'^^^^^ is a Brownian motion with drift /i + 1 started from and in 
the last step Remark 12.21 is applied. Using the scaling property of Bessel 
processes we obtain 

inf{t : R^-"^^^ = e"^} = inf{t : = 1} 

= inf{e-2^t : e'' R^J^] = 1} 

inf{t : Rt"-^ = 1}, 

where the Bessel process R!-'"^^^ is started from e^. Consequently, 

rexp(-2i?(^'))l^^,,)^^^ci.. (69) 

J 
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where Rq is distributed as with ^ exponentially distributed with pa- 
rameter 2/i. Elementary computations show that 

P{R^-^^^ >z)= z"^^, z>\. (70) 

It is interesting to notice that the right hand side of (f7n|l when extended to 
a measure on the whole of R+ can be viewed as the speed measure of 
(see, e.g., j^j A1.21 p. 133). 



5 Appendix on Ray-Knight theorems 

For an easy reference, we recall here the Ray-Knight theorems used in this 
paper (see Yor jUj and |3] for more complete statements). 

Theorem 5.1 Let B be a standard Brownian motion started from and 
local time (with respect to the Lebesgue measure) at level 

y < 1 up to Hi{B). Then the local time process {Lj^^(i?) : y > 0} is a 
solution of the SDE 



Xy = 2 r ^,df3s + 2{y M 
Jo 



in other words, 



(i) {L]j^J^^{B) : < y < 1} is a 2-dimensional squared Bessel process 
starting from 0, 

(a) {L]:^^^j(i?) : y > 1} is a 0-dimensional squared Bessel process with 
the starting value LP^^{B) obtained from (i). 

Theorem 5.2 Assume that B^^^ = and let L^(5(^)) be the local time of 
B^^^ at level y (with respect to the Lebesgue measure) up to time t. Define 
the total local time of B^^^ at level y via 

Ll^{B^^^) := lim L^(5(^)). 

t— >oo 

Then 

{L-J{B^^^) : y > 0} {Zf : y > 0}, (71) 
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and 

: y>0} {Zf : y>0}, (72) 
where Z^^'^/^), 5 = 0, 2, are solutions of the SDE 

dXt = 2^tdBt + {5- 2nXt) dt, 

respectively, with the initial value Xq exponentially distributed with param- 
eter n- In fact, the identities 171]) and I72j) hold jointly, with Z^'"^^^ = Z^^''^^'' 
but otherwise Z^^''^^^ and Z*^^'^'') are independent. 

Our final Ray-Knight theorem is for Bessel processes. The first part is 
formulated only for 3-dimensional Bessel process, and in the second part we 
take the dimension parameter 6 > 2. Let L^(i?('^)) denote the total local time 
at y of the Bessel process R^^^ (taken with respect to the Lebesgue measure). 

Theorem 5.3 a: Assume that R^^^ is started at 0. Then 

{L^^(i?(3)) : y>0} {Zf : y > 0}, 

where Z^"^^ denotes the squared Bessel process of dimension 2, started from 
0, i.e., satisfies the SDE 

dXy = 2^/X^dBy + 2dy. 
b: Assume that 6 > 2 and R^q^ = 0. Then 

where Z^"^^ denotes the 2-dimensional squared Bessel bridge (from to and 
of length 1 ). 
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